AFRL-RX-WP-TP-201 1-4218 


SITE  OCCUPANCY  OF  CHROME  IN  THE  GAMMA  PRIME 
PHASE  Ni3Al:  AN  INTEGRATED  COMPUTATIONAL  AND 
EXPERIMENTAL  STUDY  (Preprint) 


G.B.  Viswanathan  and  J.  Tiley 
Metals  Branch 

Metals,  Ceramics  &  NDE  Division 


M.  Chaudhari,  A.  Singh,  P.  Gopal,  S.  Nag,  R.  Banerjee,  and  J.  Du 
University  of  North  Texas 

Department  of  Materials  Science  and  Engineering  and  Center  for  Advanced 
Research  and  Technology 


JULY  2011 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
MATERIALS  AND  MANUFACTURING  DIRECTORATE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7750 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

Unclassified  Unclassified  Unclassified 

17.  LIMITATION 

OF  ABSTRACT: 

SAR 

18.  NUMBER  OF 
PAGES 

16 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 
Jaimie  Tiley 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

N/A 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 

Approved  for  public  release;  distribution  unlimited. 


Site  occupancy  of  chrome  in  the  gamma  prime  phase  Ni3Al:  an 
integrated  computational  and  experimental  study 

M.  Chaudhari1,  A.  Singh1,  P.  Gopal1,  S.  Nag1,  G.B.  Viswanathan2,  J.  Tiley2,  R.  Banerjee1,  and  J.  Du1 

'Department  of  Materials  Science  and  Engineering  and  Center  for  Advanced  Research  and  Technology, 
University  of  North  Texas,  Denton,  Texas  76203,  U.  S.  A. 

2Materials  and  Manufacturing  Directorate,  Air  Force  Research  Laboratory,  Dayton,  Ohio  45309,  U.  S.  A. 


Abstract: 

The  mechanical  properties  of  Ni-base  superalloys  are  strongly  influenced  by  ordered  y?-Ni3Al 
precipitates,  whose  growth  and  coarsening  are  controlled  by  the  site  occupancy  of  elements  within  the 
structures.  The  site  occupancy  behavior  of  Cr  in  y?-Ni3Al  is  investigated  by  ab  initio  based  Density 
Functional  Theory  based  computational  studies  and  3D  atom  probe  tomography.  Three  formalisms  related 
to  site  occupancy  are  discussed,  including  standard  defect  formation  formalism,  antisite  based  formalism, 
and  vacancy  based  formalism.  The  calculated  substitutional  energies  of  the  Cr  atom  at  the  Ni  and  A1 
sublattice  sites  indicate  that  Cr  has  a  strong  preference  for  the  A1  site.  Comparisons  between  the 
formalisms  indicate  that  standard  defect  formation  formalism  gives  an  inconsistent  result  and  vacancy 
based  mechanism  will  dominate  in  the  substitution  process.  The  effect  of  additional  Cr  atoms  calculated 
by  the  interaction  energies  of  two  Cr  atoms  suggests  that  Cr  atoms  prefer  to  be  on  nearest  neighbor  lattice 
sites:  either  in  the  Al,  Ni,  or  mixed  sublattice  sites.  Composition  profiles  calculated  through  atom  probe 
tomography  were  used  for  analyzing  the  site  occupancy  of  different  elements.  Based  on  the  composition 
profiles  for  Cr,  Ti  and  Al,  it  can  be  seen  that  Cr  occupies  the  Al  sublattice  site  in  the  complex  Rene-88 
alloy  and  clearly  agrees  with  the  computational  results. 

Ni-base  superalloys  find  major  applications  in  the  hot  sections  of  aerospace  engines,  industrial 
gas  turbines  and  marine  turbines  where  their  creep,  oxidation,  and  strength  properties  make  them  ideal 
materials.  The  high  temperature  mechanical  behavior  of  the  superalloys  is  critical  to  increase  the  engine 
operating  temperature  and  improve  fuel  efficiency  while  maintaining  power  generation  requirements. 

The  structure  of  most  precipitation  strengthened  Ni-base  superalloys  consists  of  a  gamma  matrix 
(y)  with  precipitates  (y’).  The  y-phase  is  a  solid  solution  with  a  face-centered  crystal  (fee)  lattice  and 
randomly  distributed  different  species  of  atoms.  The  primary  strengthening  phase  in  Ni-based  superalloys 
is  Ni3(Al,Ti),  and  is  called  gamma  prime  (y?).  It  is  a  coherently  precipitating  phase  with  an  ordered  Ll2 
(fee)  crystal  structure.  The  high  temperature  strength,  creep  and  corrosion  resistance  can  be  improved  by 
introducing  alloying  elements  such  as  Cr,  Co  and  other  elements.1,2  Partitioning  of  these  alloying 
elements  into  the  y  matrix  and  the  y’  phase  plays  a  critical  role  in  the  strengthening  mechanisms  of  these 
super  alloys. 

A  number  of  studies  have  investigated  the  partitioning  of  transition  metal  elements  (such  as  Cr, 
Co,  W,  and  Ta)  using  both  experimental  and  computational  methods.  Previous  atom-probe  tomography 
(APT)3,  scanning  electron  microscopy  (SEM)3,  atom  location  by  channeling  enhanced  microanalysis 
(ALCHEMI)  technique4,  and  other  techniques  have  concluded  that  Cr  occupies  the  Al  sublattice  in  y?. 
Atom  probe  field  ion  microscopy  (APFIM)  has  been  used  by  several  researchers5,6  to  measure  the  site 
preference  of  atoms  like  Co  in  Ni3Al.  Additionally,  several  computational  techniques  have  been  used  in 
the  literature  including  first  ab  initio  based3,7,  first  principles  method8,9,  EAM  potential  based10,11,  cluster 
variation  methods24  to  calculate  the  site  preference  for  different  alloying  additions  in  Ni3Al. 

A  lot  of  earlier  theoretical  work  to  predict  the  site  preference  was  based  on  the  standard  defect 
formalism3  which  involved  calculating  the  difference  in  formation  enthalpies  of  two  alloys,  each 
containing  one  defect  atom  on  the  Ni  and  Al  sublattice  site  respectively.  In  this  approach,  the  energy 
depends  on  the  choice  of  the  reference  state  and  is  therefore  not  complete.  Another  approach  used  by 
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Jiang  et  al 7  used  the  antisites  as  mediators  for  site  substitution.  In  this  paper,  we  present  our  results  of  a 
systematic  comparison  of  characterizing  the  site  preference  of  the  above  mentioned  formalisms  and  a 
vacancy  based  formalism  for  Cr,  which  is  the  most  common  alloying  elements  Ni3Al.  The  computational 
results  are  compared  to  3D  atom  probe  experimental  methods.  In  addition,  the  effects  of  two  Cr  atoms  in 
the  Ni-Ni,  Al-Al  and  Ni-Al  sublattice  sites  as  a  function  of  separation  distance  are  also  investigated.  All 
the  calculations  were  done  using  the  plane-wave  pseudopotential  method  as  implemented  in  Vienna  ab 
initio  Simulation  Package  (VASP). 12,13  A  plane  wave  basis  set  with  a  kinetic  energy  cutoff  of  400  eV  and 
the  Projected  Augmented  Wave  (PAW)  potentials  were  used  in  the  calculations.  The  generalized  gradient 
approximation  (GGA)  with  the  PBE14  form  was  used  for  the  exchange  and  correlation  functions.  All  of 
the  calculations  performed  were  spin-polarized.  The  atomic  structure  was  fully  relaxed  until  the  forces 
acting  on  each  of  the  atoms  were  less  than  0.01  eV/A. 

A  2x2x2  and  3x3x3  periodic  supercell  was  employed  for  defect  calculations.  A  3x3x3  supercell 
was  required  to  calculate  the  interaction  between  two  Cr  atoms.  A  k-mesh  size  of  9x9x9  and  4x4x4  was 
used  for  the  two  supercells  respectively.  The  lattice  parameter  was  calculated  to  be  3.567  Angstroms  in 
good  agreement  with  previous  reports.3 

The  formation  energies  of  intrinsic  defects  was  first  determined  in  Ni3Al  using  the  standard  defect 
formation  formalism.3  These  included  vacancies  i.e.,  VaNi  and  VaAi,  and  anti-sites  i.e.,  NiAi  and  AlNi.  The 
total  defects  involved  thermal  defect  complexes  also,  including  exchange-type  (0  ->  AlNi  +  NiAi  )  and 
Schottky-type  ( 0  3VaNi  +  VaA/)  defects.  The  formation  energies  for  the  intrinsic  defects  are  calculated 
in  both  the  supercell  sizes.  The  results  are  listed  in  Table  1,  and  are  in  good  agreement  with  the  existing 
results  of  DFT  calculations,  ’  ’  ’  experimental  studies  ’  ’  and  those  calculated  from  EAM  potentials.  ’ 

The  total  energies  of  (Ni,Cr)3Al  and  Ni3(Cr,Al)  alloy  structures  were  calculated  by  substituting 
one  Cr  atom  at  the  Ni  or  the  Al  sites.  Three  different  formalisms  were  employed  to  characterize  the  site- 
preference  behavior  of  Cr  in  y’,  including  the  standard  defect  formalism,3  antisite  based  formalism7,  and 
the  vacancy  based  formalism.  Although  the  former  two  have  been  widely  used  in  the  literature  to 
determine  site  preferences,  the  vacancy  based  mechanism  has  not  been. 

In  the  standard  defect  formalism,3  the  substitutional  formation  energies  are  calculated  as  per  the 
definition  below: 

ECr  = 

CrM 

ECr  = 

C rAl 

The  configuration  with  the  lowest  formation  energy  is  the  preferred  sublattice  for  the  alloying  atom.  This 
approach  is  not  complete  since  the  energies  depend  on  the  choice  of  reference  states  for  the  pure 
elements.  The  elemental  chemical  potentials  of  Ni,  Al  and  Cr  were  found  out  to  be  -5.748,  -3.733  and  - 
9.594  eV  atom'1  respectively.  The  energies  calculated  using  the  equation  above  in  both  the  32-atom  and 
108-atom  supercell  indicate  that  Cr  prefers  the  Ni-sublattice,  in  agreement  with  Jiang  et  al 16  and  in 
disagreement  with  Seidman.3 

The  anti-site  based  substitutional  formalism7  has  been  adopted  by  several  researchers  from  Ruban 
and  Skriver1  to  Chao  Jiang.1,7,16  In  this  formalism,  the  mediator  for  the  site  substitution  is  anti-sites.  The 

parameter  E^AI  is  defined  as  the  energy  required  in  moving  an  atom  from  one  sublattice  site  to  the 

other  sublattice  site  via  the  reaction  so  that  the  absolute  value  of  the  parameter  is  totally  independent  of 
the  elemental  reference  states  or  its  chemical  activities.  The  energy  required  in  transferring  a  Cr  atom 
from  a  Ni  sublattice  site  to  an  Al  sublattice  site  was  calculated  using  the  following  equation: 

Ni(X_^AlYCr  +  NixAlY  — »  NixAl^Y_^Cr  +  Ni^x_^AlYAl 

=  E{NixAl{Y_lfr)  +  E(Ni[x_x)AlYAl )  -  E(Ni{x_x)AlYCr )  -  E(NixAlY) 


\E{m(x_x)Cr)AlY  +  Em  )_  \ENixAlY  +  ECr  ). 
[  jjtotal  ,  |  |  jjtotal  ,  | 

\~Niy  (Ahy.\Cr)  +  EaI  /  \^NirAly  +  Ecr  /. 
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The  formation  energy  of  an  exchange  antisite  defect  in  the  Ni3Al  structure  i.e.  NiAi  +  Alm  was 
calculated  to  be  1.038  and  1.118  for  2x2x2  and  3x3x3  supercells  respectively.  If  the  value  of  E^A1  is 

less  than  zero,  then  the  reaction  prefers  going  forward,  that  is  Cr  prefers  to  go  to  A1  site,  whereas  if  the 
value  is  greater  than  the  exchange  antisite  energy,  then  the  Cr  prefers  to  go  to  Ni  site.  If  the  energy  value 
is  in  between  the  two,  then  the  Cr  atom  has  a  compositionally  dependent  site  preference. Our  calculations 
suggests  that  Cr  strongly  prefers  to  go  to  an  A1  site  in  both  cases  of  2x2x2  and  3x3x3  supercells.  This 
confirms  with  the  study  of  Jiang16  yet  shows  a  discrepancy  with  the  numbers  back-calculated  from  data 
provided  by  Seidman.3 

The  vacancy  based  substitutional  formalism  which  involves  the  mediation  of  the  vacancies  for 
site  substitution.  The  parameter  E^AI  is  defined  similar  to  the  way  it  is  defined  in  the  anti-site  based 

formalism  so  that  the  absolute  value  of  the  parameter  is  totally  independent  of  the  elemental  reference 
states  or  its  chemical  activities.  In  our  case,  we  have  studied  the  energy  required  in  transferring  a  Cr  atom 
from  the  Ni  sublattice  site  to  an  A1  sublattice  site  using  the  equation: 

Ni(X-\)AlYCr  +  NixAl(Y_i)  — ^  NixAl^Y_^Cr  +  Ni^x_^AlY 

1  —  NixAl(y_^Cr  +  Ni(X_^AlY  —  Ni^x  _^AlYCr  —  NixAl^Y_^ 

The  formation  energy  of  schottky  defects  in  the  Ni3Al  structure,  i.e.  VaAi  +  3  VaNi ,  was  calculated 
to  be  8.323  and  9.764  for  2x2x2  and  3x3x3  supercells  respectively.  If  the  value  of  E^Al  is  less  than 

zero,  then  the  reaction  obviously  prefers  going  forward,  i.e.  Cr  prefers  to  go  to  an  A1  site,  whereas  if  the 
value  is  greater  than  the  schottky  defect  energy,  then  the  Cr  prefers  to  go  to  a  Ni  site.  As  before,  the 
calculations  suggest  that  Cr  very  strongly  prefers  to  go  to  an  A1  site. 

The  calculations  show  a  discrepancy  in  the  standard  defect  formation  energy  formalism,  possibly 
from  one  of  two  sources.  One  source  can  be  the  insufficient  definition  of  formation  energy  of  the  defect  in 
the  standard  defect  formation  formalism.  The  choice  of  the  reference  state  and  the  elemental  chemical 
activity  alone  cannot  compare  with  the  size  of  the  supercell  and  may  provide  incorrect  information  on  the 
energies.  The  other  reason  originates  from  the  direct  comparison  of  the  calculated  formation  energies  of 
different  structures.  The  composition  of  the  two  structures  (yand  f)  are  different,  i.e.  both  structures 
have  a  different  total  number  and  type  of  atoms  in  them.  The  formalism  should  employ  the  same  number 
of  atoms  on  each  side  of  the  precipitate  interface  to  compare  the  formation  energies. 

Both  vacancy  and  anti-site  based  formalisms  show  that  Cr  prefers  the  A1  sublattice  site  in  f  as 
can  be  seen  from  the  negative  formation  energies.  The  vacancy  based  mechanism  has  a  larger  negative 
value  of  -1.58  eV  compared  to  around  -0.85  eV  for  the  antisite  based  mechanism,  which  suggests  that  the 
vacancy  based  mechanism  will  dominate  in  the  substitution  process.  However,  to  determine  the  total 
reaction  rate,  we  also  need  to  consider  the  concentration,  or  the  availability,  of  the  defects.  From  the 
intrinsic  defect  formation  energy  shown  in  Table  1,  the  antisite  defect  formation  energies  are  lower  than 
the  Schottky  formation  energies,  which  indicate  that  the  antisite  will  have  a  much  higher  concentration  of 
vacancy  defects.  Therefore,  both  vacancy  and  antisite  based  substitutional  mechanism  are  possible  to 
incorporate  Cr  in  the  A1  sublattice  site. 

The  effect  of  increasing  the  Cr  content  on  the  site  preference  in  f  was  investigated  by 
investigating  the  interaction  between  two  Cr  atoms  using  a  3x3x3  supercell  in  the  calculations.  The  total 
energies  of  two  Cr  atoms  were  calculated  as  a  function  of  distance  of  separation  and  the  site  occupation, 
i.e.  Ni-Ni,  Al-Al  and  Ni-Al.  A  plot  of  the  total  energies  vs  distance  of  separation  between  the  Cr  atoms  is 
given  in  Fig.  1. 

The  graph  shows  a  decrease  of  the  total  energy  as  the  distance  between  the  Cr  atoms  decreases  in 
all  three  sites  combinations,  suggesting  that  Cr  atoms  prefer  to  be  close  to  each  other  either  at  the  Ni 
sublattice  sites,  the  A1  sublattice  sites  or  in  the  mixed  sublattice. 
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The  results  of  the  DFT-based  calculations  of  Cr  site  occupancy  in  the  y?-Ni3Al  phase,  exhibiting 
the  ordered  Ll2  structure,  in  a  ternary  Ni-Al-Cr  alloy  have  been  experimentally  validated  in  a  commercial 
nickel  base  superalloy,  Rene  88  (55.9Ni-17.9Cr-4.5Al-12.8Co-4.5Ti-2.4Mo-1.2W-0.4Nb,  all  at%).  Rene 
88  is  a  complex  alloy  with  at  least  five  major  alloying  elements,  Al,  Ti,  Cr,  and  Co,  amongst  which  A1 
and  Ti  primarily  partition  to  the  ordered  y?  phase  while  Cr  and  Co  partition  to  the  disordered  y  matrix.17,18 
Since  the  composition  of  Cr  is  the  highest  amongst  all  these  alloying  elements  in  the  Rene  88  alloy, 
coupled  with  the  fact  that  there  is  an  approximate  3at%  Cr  composition  within  the  y?  phase,  the  Cr  site 
occupancy  is  of  substantial  interest  and  has  been  investigated  using  atom  probe  tomography  (APT).  The 
alloy  sample  used  for  this  study  was  solutionized  at  1150°C  for  30  mins,  continuously  cooled  to  room 
temperature  at  a  rate  of  24°C/min,  and  subsequently  aged  for  50  hours  at  750°C.  The  inset  in  Fig.  2(a) 
shows  an  energy-filtered  TEM  micrograph  exhibiting  a  bimodal  size  distribution  of  y?  precipitates  (darker 
contrast).19  Samples  for  APT  investigations  were  prepared  using  standard  eletro-polishing  techniques  as 
described  elsewhere.20  Fig.  2(a)  also  shows  a  section  of  the  APT  reconstruction  (only  Al  atoms  plotted) 
with  the  y?  regions  exhibiting  a  higher  Al  concentration.  Within  this  three-dimensional  (3D) 
reconstruction,  a  primary  y?  precipitate  (marked  in  Fig.  2(a)),  whose  <001>  direction  lies  along  the  Z-axis 
(axis  of  reconstruction),  was  selected  for  analysis  since  the  highest  spatial  resolution  is  afforded  along  this 
direction.20,21  The  (001)  ordered  planes  of  Al,  characteristic  of  the  ordered  Ll2  structure,  are  clearly 
visible  within  the  y’  precipitate  in  this  reconstruction.  It  should  be  noted  that  these  planes  are  not  perfectly 
aligned  and  suffer  from  errors  associated  with  the  reconstruction,  possibly  leading  to  artifacts  such  as  the 
presence  of  Al  atoms  between  the  ordered  (001)  planes,  not  allowed  in  the  y?  crystal  structure. 
Additionally,  the  presence  of  Al  antisite  defects  on  Ni  sites  cannot  be  ruled  out.  Composition  profiles, 
averaged  across  a  cylinder  passing  through  the  y?  precipitate  (shown  in  Fig.  2(a)),  perpendicular  to  the 
(001)  Al  planes,  were  used  for  analyzing  the  site  occupancies  of  the  various  elements.  The  composition 
profiles  for  Ni  and  Al  exhibit  periodic  modulations  as  shown  in  Fig.  2(b)  where  the  peaks  of  Ni  align  with 
the  troughs  of  Al  and  vice  versa,  validating  the  presence  of  alternating  Ni  and  50%Ni+50%Al  sublattices 
sites  along  the  <001>  direction  of  y\  Similar  composition  profiles  for  Al  and  Ti  have  been  plotted  in  Fig. 
2(c).  Previous  reports  in  the  literature16,22  clearly  indicate  that  Al  and  Ti  occupy  the  same  sublattice  site  in 
the  ordered  Ll2  structure.  While  the  profiles  for  Al  and  Ti,  shown  in  Fig.  2(c),  indicate  that  these  two 
atomic  species  nominally  lie  on  the  same  sublattice,  there  is  a  clear  marginal  shift  of  the  Ti  peaks  to  the 
left  of  the  corresponding  Al  peaks.  This  shift  can  be  attributed  to  a  reconstruction  artifact  arising  from  the 
difference  in  the  evaporation  fields  and  time  of  flight  (TOF)  associated  with  these  two  elements.20  Thus, 
while  the  evaporation  field  for  Al  is  -  19  Vnm'1,  the  value  for  Ti  is  substantially  higher  (~  41  Vnm'1). 
Finally,  based  on  the  composition  profiles  for  Ti  and  Cr,  shown  in  Fig.  2(d),  it  can  be  concluded  that  Cr 
occupies  the  same  sublattice  site  as  that  of  Ti  (and  Al)  in  the  ordered  Ll2  structure.  Therefore,  even 
though  Rene  88  is  a  significantly  more  complex  alloy  compared  to  the  ternary  Ni-Al-Cr  system,  Cr 
occupies  the  Al  sublattice  in  this  commercial  alloy,  in  agreement  with  the  results  of  DFT  based 
calculations  presented  in  this  letter. 

To  summarize,  the  Cr  atom  shows  a  strong  preference  to  the  Al  sublattice  site  from  the 
composition  profiles  of  atom-probe  data  and  ab  initio  based  calculations.  The  defect  formation  formalism 
shows  inconsistent  results  because  of  the  incomplete  definition  and  choice  of  reference  states.  Both 
vacancy  and  antisite  based  formalisms  show  a  strong  preference  for  Al  site  with  the  prior  being  a 
dominant  substitutional  process.  The  interaction  between  two  Cr  atoms  shows  that  the  two  Cr  atoms  will 
tend  to  be  as  close  as  possible  to  each  other. 
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Table  1  Intrinsic  defect  formation  energies. 


Designation 

Defect  type 

Defect  formation 

energy 

2x2x2  3x3x3 

Other  Studies 

VaA1 

A1  vacancy 

3.647 

3.570 

3.097,  2.658 

1 

> 

Ni  vacancy 

1.559 

2.064 

1.157,  1.87s,  1.609,1.8023 

Mai 

Ni  antisite 

2.096 

1.957 

0.9863,  2.047 

AlNi 

A1  antisite 

-1.058 

-0.838 

0.7423,  -0.927 

0  ->  AlNi  +  NiAi 

Exchange 

1.038 

1.118 

1.7293,  1.127,  1.448,  1.029,  1.1510,  1.6711 

0  3VaNi  +  VaAi 

Schottky 

8.323 

9.764 

6.547,  8.268,  6.509,  6.3310,  6.7311 

23  -  Experimental  study 


3,16  -  ab  initio  DFT  method 
8,9  -  First  principles  method 
10,11  -  EAM  potential  method 


Table  2  Site  preference  energies  calculated  from  the  standard  defect  formalism,  antisite  and  vacancy 
mechanism  with  different  super  cell  sizes. 

Standard  defect  formation 

Our  Calculations 

Other  Studies 

formalism  -Cr 

2x2x2 

3x3x3 

Seidman3 

Jiang16 

77 Al 
hCr 

1.193 

1.363 

0.565 

1.33 

77  Ni 
hCr 

0.946 

1.212 

0.648 

1.29 

Our  Calculations 

Other  Studies 

Antisite  based  formalism-Cr 

2x2x2  3x3x3 

Jiang3 

Seidman16 

ttt  Ni — }Al 
^ Cr 

-0.810  -0.917 

-0.50 

0.695 

Vacancy  based  formalism 

V  dL/dlluy  UdoCLl  lUIllldllMIl“v^I 

2x2x2 

3x3x3 

pNi^Al 
^  Cr 

-1.580 

-1.585 
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Cr-Cr  in  Ni-Ni  site 
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|  -618.8000 
5*  -619.0000 
^  -619.2000 

OX) 

<5  -619.4000 
-619.6000 


Distance  between  two  Cr  atoms 
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Cr-Cr  in  Ni-AI  site 


B  -620.7500 
|  -620.8000 
>•  -620.8500 
~  -620.9000 
2?  -620.9500 
|  -621.0000 
-621.0500 


2  3  4  5  6 

Distance  between  two  Cr  atoms 


Cr-Cr  in  AI-AI  site 


E 

Q 

-622.2 

-622.25 

-*-> 

rz 

-622.3 

<D 

-622.35 

>* 

OX) 

-622.4 

i— 

<1> 

c 

-622.45 

LU 

-622.5 

3  4  5  6  7 


Distance  between  two  Cr  atoms 


Fig.  1  Interaction  energy  as  a  function  of  Cr-Cr  distance  in  different  sublattices. 
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Fig.  2  Atom  probe  results:  a)  data  region  for  the  ATP  data  collection,  inset  shows  an  energy  filtered  TEM 
image  of  y  (lighter)  and  y’  (darker)  phases,  b)  plot  of  the  Ni  and  A1  at%  compositions  from  ATP,  c)  plot 
of  the  A1  and  Ti  at%  compositions  from  ATP,  and  d)  plot  of  Ti  and  Cr  at%  compositions  from  ATP. 
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